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Abstract: Within the framework of the Lee Wick Standard Model (LWSM) we investigate 
Higgs pair production gg — > hoho, gg — > h^pQ and top pair production gg — > it at the 
Large Hadron Collider (LHC), where the neutral particles from the Higgs sector (ho, ho 
and po) appear as possible resonant intermediate states. We investigate the signal gg — > 
hoho — > 6677 and we find that the LW Higgs, depending on its mass-range, can be seen 
not long after the LHC upgrade in 2012. More precisely this happens when the new LW 
Higgs states are below the top pair threshold. In gg — > it the LW states, due to the 
wrong-sign propagator and negative width, lead to a dip-peak structure instead of the 
usual peak-dip structure which gives a characteristic signal especially for low-lying LW 
Higgs states. We comment on the LWSM and the forward-backward asymmetry in view 
of the measurement at the TeVatron. Furthermore, we present a technique which reduces 
the hyperbolic diagonalization to standard diagonalization methods. We clarify issues of 
spurious phases in the Yukawa sector. 
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1 Introduction 



1.1 The Lee- Wick Standard Model 

The investigation of the mechanism of electroweak symmetry breaking (EWSB), responsi- 
ble for the generation of fermion and gauge boson masses, is one of the primary tasks of 
the Large Hadron Collider (LHC) at CERN. The scalar Higgs particle realizes this mecha- 
nism in the Standard Model (SM), in a rather efficient way, at the expense of divergences 
quadratic in the cut-off. The latter fact, known as the hierarchy problem, is taken as 
an indication of the incompleteness of the SM and is at the heart of many models be- 
yond the SM (BSM). An example of which is the Lee- Wick SM (LWSM) [1] where ideas 
to soften ultraviolet (UV) divergences in QED from the seventies [2, 3] were extended to 
chiral fermions and non-abelian gauge theories [1]. Most importantly it was shown that 
the LWSM is renormalizable and free from quadratic divergences [1] thus joining the list of 
models addressing the hierarchy problem successfully. In LW field theories higher derivative 
(HD) terms are added and terms quadratic in the fields are resummed into the propagator 
rather than treated as perturbations, ameliorating the UV behavior of perturbation theory. 
This results in additional poles in the propagators for which auxiliary fields (AF) can be 
introduced to cast the theory in terms of interactions with mass dimension no greater than 
four 1 . The additional fields are interpreted as LW partner states and do have the wrong- 
sign propagator, aka Pauli-Villars regulators. The key idea of Lee and Wick is that the 
LW ghost particles never appear as asymptotic states in detecters, nowadays reminiscent 
of the Faddev-Poppov ghosts in non-abelian gauge field theories. The connected issues of 
unitarity and causality which were debated in the seventies, e.g. the Erice lectures [4, 5], 
and reconsidered recently in [6]. Most notably the width becomes negative and requires a 
deformation of the contour to avoid new cuts [7] which assure no new asymptotic states. 
The status of LW field theories is that there are no known counterexamples to unitarity in 
perturbation theory up to today and that causality can be violated but only at distances 
below M^j. It has been suggested that the violation of causality can be tested at the 
LHC [8]. The usual non-perturbative formulation via the path- integral seems difficult [9] 
but recently a restrictive path-integral was proposed where the contour prescription can 
be derived [10, 11]. 

Further conceptual issues of phenomenological nature have been investigated such as 
the behaviour at high temperature [12], unitarity of massive LW vector boson scattering 
[13], the compatibility of the see-saw and the absence of quadratic divergences [14], the 
running of couplings [15], UV-properties of LW field theories [16], even higher derivative 
LW field theories [17, 18] and LW fields and gravity [19]. The cosmology of LW field theories 
has been investigated in [20]. Phenomenological studies include LHC and linear collider 
signals of LW gauge bosons [21, 22], flavour changing neutral currents [23], electroweak 
precision observables (EWPO) have been investigated in [24] and [25] where gauge boson 
and fermion masses are found to be constrained up to a few TeV. 

It is amusing that in the AF-formalism the LWSM seems fine tuned with respect to the hierarchy 
problem whereas this is not the case in the HD-formalism as a single term is added in each sector. 
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1.2 The Higgs-sector of the Lee- Wick SM 

The LW Higgs sector has been investigated in [26-29]. The neutral part consists of the 
CP-even ho, ho, which are the SM-like and the LW-like Higgs boson, and the CP-odd LW- 
like scalar po- The SM as well as the LW Higgs sectors are not easy to constrain, neither 
indirectly through loops nor directly through signals. First for large Higgs masses the 
latter enters only logarithmically, rather than quadratically, at one-loop [30]. Second the 
Higgs couples via Yukawa terms to fermions and is therefore highly suppressed in di-lepton 
signals h — > l + l~ 2 . 

A salient feature of the LWSM, at least in its minimal version [1] , is that there's roughly 
a single new parameter per sector. It's the mass in the HD formalism which predicts all 
masses and couplings in the language of the AF formalism. In this respect the LWSM 
resembles so-called sequential SM extensions. The aim of this paper is to investigate the 
effect of a low lying Higgs sector, as a function of this single new parameter and the Higgs 
mass. We focus on channels, accessible at the LHC, where the additional Higgs appear as 
intermediate states at or close to resonance. 

• Higgs boson pair production is beyond reach at the LHC in the SM [32] . In extensions 
of the SM its a different quest as particles, with appropriate couplings and masses 
above the two Higgs threshold, can enhance the cross section by orders of magnitude 
without contradicting current constraints 3 . We consider gg — > hoho and gg — > poho- 
We find that the cross section of the latter can be enhanced by roughly three orders 
of magnitude with respect to the SM for a sizable range of masses. That is to say if 
the LW Higgs is above the SM-like Higgs pair threshold and not to far above the top 
pair threshold, 2mh < m^ ^ 1.5(2m t ). If the latter bound is approached top pair 
production becomes the main channel: 

• top pair production through gluon fusion does not suffer from low cross sections and 
has already been observed at colliders. The cross section of the invariant mass of the 
top-pair Mu has been identified as an attractive observable to see resonance effects 
through interference with the QCD-part a long time ago e.g. [34]. LW field theories 
have a very different pattern in that the wrong-sign propagator and width lead to a 
dip-peak rather than a peak-dip structure in the spectrum. It should be added that 
such effects can and do also appear in strongly coupled theories such as low energy 
QCD as discussed in section 4. 

The paper is organized as follows: In section 2 we give an overview of the Higgs 
and quark sectors within the LWSM. In sections 3 and 4 we discuss Higgs pair and top 
pair production from a theory point of view. In section 4.1 we comment on the top 
forward-backward asymmetry in view of the current TeVatron results. In section 5 we 

2 This is why, in our opinion, the LW Higgs is not a candidate for the VKjj'-excess at the TeVatron [31] 
as it should already have been seen in VKZZ-signal or Wfcb-signal. 

3 A well-known example is minimal supersymmetric SM (MSSM) [33]. In fact the LWSM Higgs sector 
particle content corresponds to a type-II two Higgs-doublet model with a new LW Higgs mass scale as a 
single new parameter with tan j3 = 1. The masses of the different Higgs particles are discussed in section 2. 
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present plots. In section 6 we investigate the signal gg — > h^ho — > bb^. In section 7 
we conclude. In appendices A and B we present further details of amplitudes for Higgs 
pair and top pair production, respectively. In appendix C a method that reduces the 
hyperbolic diagonalization to standard techniques is presented. In appendices C.l and C.2 
we present tree-level mass sum rules. Further, we clarify the issue of spurious phases versus 
CP-violating phases in the fermion mass matrices. 



2 The Lee Wick Standard Model 

We shall discuss the Higgs and Yukawa sectors directly in the auxiliary field formalism and 
refer the interested reader to [1] for the connection with the higher derivative formalism. 

2.1 Higgs sector 

The Lagrangian of the Higgs sector in the auxiliary field formalism assumes the following 
form [1]: 



C = (D^H)\Df"H) - (D^H)\D^H) + M 2 H H ] H - V(H - H) , 



(2.1) 



where D„ = 8^ + i(A„ + A M ) with A M = gA^T a + g 2 W*T a + 9l B^Y for SM gauge 
fields and analogously for the LW gauge boson for A^. The Higgs potential is V(H) = 
\/4,{H^ H — v 2 /2) 2 . The mass Mh is the mass scale of the higher derivative LW mass scale. 
In the unitary gauge the two doublets are 



H T = [0, (v + h )/V2] , H T = [h+, (h + ipa)/V2] . 



(2.2) 



It is worthwhile to emphasize that, prior to mixing, the SM but not the LW CP-even 
neutral Higgs acquires a vacuum expectation value: 



(h ) = v , (h ) = . 



(2.3) 



We note the standard abuse of notation in not denoting the massless as well as the massive 
Higgs field by h$. With (2.2) the mass Lagrangian assumes the following form: 

A 2/ , r , 2 Ml 

Anass = {ho - ho) + -y 

We note the mixing between the Higgs scalar and its LW-partner. The neutral CP-even 
Higgs field can be diagonalized by a symplectic rotation: 



}v 2 (h - h ) 2 + ^(hoho+paPo + 2h+h-) 



(2.4) 



(cosh 4>h smh(f>h 
smh(f>h cosh^ 

for which the masses of the Higgs sector are given by, 




(2.5) 
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and for completeness we have indicated the CP quantum numbers as well. For obtaining 
Feynman rules in terms of the physical masses the following relations are useful [26]: 



2m Lphys 



and 



2.2 Yukawa Interactions 



s H = cosh^ 
cosh 6h 



r h 



1 



m 



(2.7) 



/i ,phys 



(1-^ ) 1/2 ' 

sinh 4> h 



(1-^ ) 1/2 ' 



(2i 



In order to discuss the Yukawa terms, it is helpful to first discuss the fermions. We shall 
closely follow ref. [26]. However, we choose a slightly different basis for the fermions and re- 
fer the reader to appendix C where a method is outlined how the hyperbolic diagonalization 
can be performed using standard tools. 

The kinetic term of the AF Lagrangian is given by: 



with 



iT 



(T L ,t' L ,T L ), ^ = (t R ,i R ,T' f 



R) ! 



(2.9) 



(2.10) 



where all capitalized components are part of an SU(2) doublet; e.g. Ql = (Tl,Bl) t . It 
is noteworthy that a chiral fermion necessitates two chiral fermions which in turn form a 
massive Dirac fermion. This becomes explicit in the basis chosen above 



Mtm 



( 



\ 



nit —m t 
-m t -M u m t 
-Mr 



(\ o o \ 

0-10 

V° -v 



(2.11) 



which differs from the one in [26] . Note though that all physical masses remain unchanged 
under change of basis. The mass matrix is diagonalized by symplectic rotations Sl and 
Sr: 



*£CR),phys = mS[ {R) V3^ L(R) , M,phys??3 = S R M t r] 3 S L 
which leave the kinetic terms invariant by virtue of 



ShmS ] L = % and S R i] 3 S R 



(2.12) 



(2.13) 



Now we may turn to the Yukawa sector for which we only write down the neutral Higgs 



part: 



ip )[^ R g t ¥ L -^ L g\¥ R ) , (2.14) 



- 5 - 



where the g matrix has non-diagonal entries which allow for transitions between LW- 
generations and is given in the initial and physical basis by: 



fjt 



' mt —nit^ 

—nit mt 
^ i 



9t,ph ys = S J R gtS L . ( 2 -15) 



3 Higgs boson pair production 

We shall parametrize the gg — > /lo^o matrix element as follows: 

M(gg -»• h ho) = 3^2^^ ( A ° P ° + e(pi)g e{p 2 )l , (3.1) 



with analogous conventions for gg — > poh^. The pre-factor arises as follows: l/25 ab due to 
the colour trace, 1/4 from perturbative expansion, the fraction g 2 jv 2 from the couplings of 
the vertices and l/(4-7r 2 ) is factored out in order to give simple results for the amplitudes. 
The parity-even projectors on gluon spin and 2, Po and P2, as well as their parity-odd 
counterparts, Po and P%, are defined in appendix A. The parton cross section for 2 — > 2 
scattering process for two massless incoming particles is given by l/(167rs 2 )|A^ | 2 [36] and 
averaging over initial state polarizations 1/4 and colour 1/(N 2 — l) 2 = 1/64 one arrives 
at 4 : 

dajgg ->■ h h ) 1 1 9 4 2 2 

= ^19^^V" A ° + ^2 J (3.2) 

This result is for identical particles. In the case the particles in the final state are not 
identical one has to multiply by a factor of two 5 . The spin amplitudes, parity-even 
and odd, receive contributions from the triangle and box diagrams, c.f. figure 1 (left) and 
(right) respectively, whereas the spin 2 amplitudes only receive contributions from the box 
diagrams: 

A = a£ + A 1 , A 2 = A% . (3.3) 

For what follows it is important to notice that the gluon-quark vertex is diagonal in 
LW-generation space whereas the Higgs-quark vertex is not (2.14). Since, the Higgs-quark 
vertex does not contribute to the triangle graph the latter can be obtained from the SM 
with simple corrections for vertices as described in appendix A.l. The modification of the 
box graphs are twofold. First, the external Higgs particles are modified by the mixing factor 



s- 2 
H 

vertex the LW-generations mix (2.14) as discussed above. We find that these modifications 



s 2 r - r as for the triangle. Second, one has to take into account that at the Higgs-quark 



4 This agrees with [37] with the following identifications: |.4o| 2 = |gaugef | 2 and \A2^ — |gauge2| 2 at the 
difference that here Ao,2 are meant to include the LW contributions as well. 

5 We have thus implicitly assumed that the variable t is understood to be integrated over its entire 
domain despite the Bose symmetry in the identical particle case. 
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(a) (b) 

Figure 1. (a) Triangle graphs for q — (t,i,T,b,b, B) and (b) one out of six box graphs for 
qi,qj = (t,t,f,b,b,B). 

are most efficiently presented as follows: 

3 

"4o (gg -> h h ) = s 2 H _ & ^ (fu(i,j) + f 55(1,3)) 

3 

^0 (gg -> ^oPo) = ^ (/ib(*> + /5i(*. J*)) . ( 3 - 4 ) 

where 

flavours 

fxx(i,j)= S(X) f ] ij [ri 3 S(X) t ] ji (a ) xx (rni,rn j )] 
f 

flavours 

fxv(i,j)= £ [%5(X) / ] ii [^P(y) / ] ii (ao)^y(mi J m i )]. (3.5) 
/ 

In regard to the formulae (A. 17) it is important to notice that the /io and po are associated 
with the the momenta p^ and p^ respectively as can be inferred from the formula in appendix 
A. 2. The couplings S(P)x,Y, which follow from eq. (2.14), are: 

S(l) t = ±(gt + g t ), S(5) t = ±(gt-g t ), 

P(l)t = \(- 9 t+9t), P(S) t = \(-gt-9t), (3.6) 

where the top flavour was chosen as a representative and the subscript phys has been omit- 
ted on the Yukawa couplings for the sake of notational brievety. The 773 = diag(l, —1,-1) 
matrices take care of the signs of the SM and LW propagators respectively and the cou- 
plings gx,Y govern the LW-generation transitions. The spin 2 structures and ^2 are 
completely analoguous. 
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4 Top pair production 



In this section we discuss the interference between the QCD background and resonant 
particles in top pair production in a qualitative manner. 6 In the LWSM potential resonant 
particles, that couple to the top triangle loop and decay into top pairs are the ho, ho, po, Z 
and Z corresponding to the diagrams shown in figure 1(a) and figure 12(a,b,c), respectively, 
with the Higgs final states replaced by top pairs. Here we shall neglect the Z and the Z 
as the former is far off-shell at s > 2mt and the latter is severely constrained by di-lepton 
searches to be heavier than 1 TeV and by electroweak precision measurement to be in the 
multi-TeV range. The corresponding amplitudes, which consist of triangle graphs only, are 
easily obtained from the one for Higgs-production and are given in appendix B. 

The interference effect of an intermediate resonance gg —> R — > tt, where R stands for 
a generic resonance, takes the following form [34]: 



~^(99 -> tt) I interference = -|c(s)|Re 



A 



m% + im R T R 



R 

= -\c(s)\ ((a - m|)Re[Z A ] + m R T R lm[l A }) , (4.1) 

where l/\ = l/\{s/Am1) is the appropriate triangle loop function, c(s) is a well-known 
function of s [34], c(s) differs from c(s) by a constant and s the invariant mass of the two 
gluons entering the loop. If there is no loop function then the term above will lead to a 
peak-dip structure passing from constructive to destructive interference at s = m R . The 
loop-function does not change this pattern in the case where the resonance is a scalar or 
a pseudoscalar [34] as the real and imaginary part of the loop function are positive. The 
pattern persists for a spin-1 particles as well as can be inferred from the plots in reference 
[44]. Thus the question what happens in the LW case. Due to the negative sign of the 
propagator and the width, 



-J-{gg ->■ tt)|LW-mterference = -|c(s)|Re 

as 



-Z A (s/4mf 



(s - m 2 R ) - im R T R 



R) 

= -\c(s)\ (-(s - m|)Re[Z A ] + m R T R lm[l A }) , (4.2) 

the (s — mJj)Re[Z/\] term flips sign . Assuming that neither the width nor the imaginary 
part of the loop function Z A are anomalously large, this leads to a dip-peak structure. In 
fact the passage from destructive to constructive interference, which we shall call M R , does 
not coincide with the exact location of the resonance: 

Examples of the effect are shown in figure 2. The dip-peak structure is a unique feature 



6 We note that in ref. [35] the authors explored these types of interferences in the context of minimal 
supersymmetric standard model and Little Higgs models. 

7 It is crucial that the intermediate resonance couples to the tops from the loops and the final state tops 
as otherwise a minus could be absorbed in either one of the couplings. 
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Figure 2. The cross section <r(gg — > it) as taken from [34] qft with energy dependent width. The 
solid line is the LO QCD contribution. The dashed(dotted) lines correspond to a resonance mass 
m,R = 400(600) , GeV. The left(right) hgure corresponds to the usual (LW) resonance-type. 

of LW field theories, produced via gluon fusion through the top triangle, in the case of 
a well isolated resonance. We would like to add to that in the case where the masses of 
two resonances are close to each other their mixing has to be taken into account by the 
so-called K-matrix formalism e.g. [45]. It is important to realize that a dip-peak structure 
is present in the 7r-7r-scattering spectrum for the /o(980) meson due to the extremely broad 
/o(600) (cr-meson) [46]. Thus strongly coupled extensions of the SM, such as technicolor, 
might have similar signals as the LWSM. 

4.1 A comment on the top forward backward asymmetry 

Currently, the top forward-backward asymmetry (tAFB), A^ B = 0.475(114) for M# > 
450 GeV at [47] at 5.3 ftT 1 , deviates from SM prediction A l j B = 0.088(13) [47] at about the 
3<7-level at the TeVatron 8 . The SM prediction originates from a charge asymmetry which, 
due to the fact that the TeVatron is a pp-collider, translates into a forward-backward 
asymmetry. Thus the question is whether the LWSM has the potential to explain this 
discrepancy. A nice summary of perturbative approaches to the tAFB is given in reference 
[49]. The LWSM qualifies at the same level as a Z'-model with SM-like couplings, where 
the role of the Z' is taken by Z. 9 We roughly get Ap B ~ 0.01, for a = 1 TeV, at best 
which is in the right direction but too small to join into the current excitement. 10 Note, 
as only the absolute value of the propagator enters, the wrong-signs of the propagator and 
the width do not matter. We have not evaluated the interference of the Z and the SM Z, 
but expect it to be of similar size. 

5 Numerical results 

We compute cross sections for pp — > /lo/io/Po^o an d the differential cross section pp —> it via 
gluon fusion at the LHC for V~S = 7/14 TeV respectively 11 . We denote the pp center of mass 

8 The very recent DO-results at 5.4 fb -1 is much closer to the SM value [48]. 

9 The LWSM does not qualify as an axi-gluon, nor are there large flavour changing couplings between 
the first and third generation in its minimal version. 

10 Note 171% = 1 TeV is even a bit low in regard to electroweak precision data [24]. 

11 The cross section for vector boson fusion qq — > hohojj is about 2% of gg — > hoho and thus negligible. 
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energy by capital S and the partonic center of mass energy by lower case s throughout this 
paper. The renormalization and factorization scale has been chosen to be \i T = \x$ = 2m^ 
for pp — > hoho and /x r = fif = m/j + mp for pp — > poho- We use the MSTW 2008 LO 
(90% C.L.) for parton distribution functions with the strong coupling calculated to one- 
loop order for a s (mz) = 0.13939 [50]. We use LO predictions for gg — > hoho/poho. The 
NLO corrections in the later case are rather large [51]; almost 100% as can be inferred from 
figure 6 of that reference. Fortunately, the shape of the corrections are almost identical to 
the LO result and thus should not distort the analyses too much. For gg — > it we also use 
LO predictions with the factorization scale set to \it = mt and the renormalization scale 
set to \i r = mt(m^) for gqq (gg<p for <f> = po,ho,ho) couplings for which we comment in 
section 4. 

For the numerical computations we have used various computer packages to be re- 
ferred to below. The FeynArts [52] model file has been generated automatically using 
LanHEP [53]. The resulting model files were modified to allow for wrong-sign propagators 
in the auxiliary field formalism. Fortran code for the cross sections was generated with the 
use of FormCalc [54] . All loop integrals were computed using LoopTools [54] . 

The width-mass ratios, widths and branching ratios for ho and ho are depicted in 
appendix A. 5. in figures 13,14(left) and 15 respectively. They will be referred throughout 
and serve to understand the results qualitatively. Possibly the most important aspect 
for further understanding is that the width of the ho (in figure 14(left) appendix A. 5) 
raises significantly when the it-threshold is crossed (in parameter-space > 2mt) and 
is relevant for the triangle diagrams with intermediate ho. 

5.1 Contraints on LW mass scales 

Before presenting the main results the new LW scales have to be discussed. There are six 
LW mass scales plus the mass of the SM-like Higgs boson out of which five are constrained 
to be rather high and generally do not impact on our investigation The parameters are: 

• The scales M\ and M2 of the LW gauge bosons associated with U(l)y and SU(2)x 
are constrained by electroweak precision measurements to be in the multi-TeV range 
[24] . We shall set Mi = M2 = 1 TeV throughout this paper as in this range the 
masses have no major influence on our results. 

• The fermion mass scales Mq, M u and 12 are constrained through loop-contributions 
to electroweak precision measurements to be in the multi-TeV range [25]. For the 
gg — > hoho channel the fermion mass scale has little influence for Mq = M u = Mj > 
500 GeV as can be inferred from the appendix A. 5 figure 16. There are no qualitative 
changes when one goes away from the limit of equal masses and we therefore assume 
the the fermion mass scales to be 500 GeV in the plots. For the gg — > poho channel 
there are some threshold effects due to the box diagrams. 

12 Due to chiral suppression, only heavy flavours are relevant. This statement can be inferred from the 
HD formalism. Thus only the top and the beauty quarks are taken into consideration. 
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• The masses of the two neutral CP-even Higgses ho and h,Q (2.6). Every other param- 
eter, in the Higgs sector, can be expressed in terms of these two. In particular the 
pseudoscalar mass, at tree-level satisfies (2.6) 

m Po = m h + m l . (5- 1 ) 
where we have dropped the subscript "phys" and shall do so in the remainder of 
this paper. The Higgs parameter-space has already been studied in other works. 
The collider analysis of gg —> ho — >• 77 [26] was extended to final state channels "fZ 
and WW in [29]. A part of the parameter-space has been found to be excluded by 
TeVatron results, c.f. figure 3 in that paper. It has to be added that this work was 
done in the narrow width approximation. Inspecting the plots in figure 13 it would 
seem that the effect of the width should be moderate in most of the parameter-space 
that has been excluded. The overlap of the interesting parameter-space and their 
excluded region is rather small and we leave it to the reader to convince him or 
herself of this fact. Using the correspondence of the the LWSM Higgs-sector and the 
type-II two Higgs doublet model mentioned, in the introduction, the effects of the 
charged Higgs boson h + on flavour physics were investigated in reference [27]. Using 
NLO predictions for b —> S7, neglecting the influence of all other LW states, which is 
consistent with our analysis, it was found that > 463 GeV at the 95% confidence 
level. Together with the tree-level relation (5.1) and mp = m,h + (2.6) this sets a 
significant constraint on the lower range of our parameter space. Concerning this 
indirect bound there are two remarks to be made. First, the individual theoretical 
uncertainties were added in quadrature, which is common practice, and thus the 
uncertainty might be considered to be a little bit on the low side. Second, the tree- 
level relation between the Higgs masses might receive significant radiative corrections 
due to the large top mass which is the case in the MSSM. 

We would like to add that the limit of degenerate masses of the ho and ho, parametrized 
by r/j = TUho/m^, is somewhat delicate. In connection with real particles, in the sense 
of parton level, it does not make sense to treat them separately. This can be seen in the 
pole in rh in s H _jj = (1 + rj^)(l — r^) -1 / 2 . For virtual particles it is best to resort to 
the HD-formalism where everything should remain consistent. In regard to these points 
we disregard the parameter space where 

r ho > 0.8 , r ho = ^ , (5.2) 
m ho 

which is somewhat more conservative than the value rh Q > 0.9 chosen in [8]. It would be 
interesting to study these effects, from scratch, in the HD-formalism and find the relation 
to the K-matrix formalism [45] used to improve on two nearby Breit-Wigner resonances in 
usual field theory. 

5.2 Results for gg — > /io^o 

The main point is that for slightly above 2m/i the cross section is large, three orders 
of magnitude larger than the one of the SM, dominated by the resonant contribution in 
the triangle graph figure l(left). This is reminiscent of the situation in the MSSM [33]. 
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Figure 3. The cross section (in fb) for gg — ¥ hoho via gluon fusion at the LHC for y/S = 7/14 TeV 
respectively versus the mass of the h , , for three different values of to/, . 



7 TeV 



14 TeV 



600 
500 
400 
300 
200 
100 




Log[tr(fb)] 

I 



600 

500 ; 
400 ; 
300 ; 
200 ; 

100 U 





200 300 400 500 600 700 800 



200 300 400 500 600 700 800 

mr 



Figure 4. Contour plot of the total cross section (if fb) for gg — > hoho (y/S — 7/14 TeV respectively) 
versus the light Higgs boson mass, mh and heavy Higgs boson mass, for Mq = 0.5 TeV and 
Mi = Mi = 1 TeV. Note figure 3 corresponds to horizontal sections in this plot. 



In figure 3 we show the total cross section for gg — > hoho for y/S = 7/14 TeV respec- 
tively as a function of for three different values of m^ . More detailed information can 
be inferred from the contour plots in the (mh , )-plane shown in figure 4. As mentioned 
above one observes a sharp raise of the cross section when the LW Higgs mass crosses the 
threshold 2mh , c.f. figure 3. For higher the resonance contribution decouples and 
finally approaches the SM value. An interesting effect arises when the top threshold is 
reached. For the observation to be made below recall that the process is dominated by the 
triangle graph with an intermediate LW Higgs propagator of the form (s— m~ — im^T^) . 
The slight dip in the branching ratio, c.f. figure 15(right), below the W-threshold results 
in a slight raise of the curve in case the where m^ < rnt- Once the it-threshold is reached 
the rapidly growing decay rate is damped through the additional relevant part in the prop- 
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Figure 5. The cross section (in fb) for pp —> hopo via gluon fusion at the LHC for VS = 7/14 
TeV respectively versus the mass of the ho, m h , f° r three different values of m^ . Note the 
kinks are due to crossing thesholds in corners of the box diagrams as described in the text; recall 
M Q = M u = M d = 500 GeV. 

agator. Note the lower part of the blue curve raises. In the HD-formalism this can be 
understood by the to the two poles rrih and approaching each other. 

5.3 Results for gg — > hofto 

The cross section for y/S = 7/14 TeV respectively with fixed are shown in figure 5. 
The corresponding contour plots are shown in figure 6. The crucial difference to gg — > 
ho — > hoho, in terms of the triangle diagram, is that there's no parameter region where 
there's a dominant resonance effect. The diagrams are shown in figure 12: the intermediate 
Z and Z are either too light or too heavy respectively and the process gg — >• po — > poho 
is not on-shell. There's a remnant of the latter effect when the mh is relatively small 
and po ~~ > Poho approaches an on-shell configuration. The cross section is enhanced for 
r h — > 0.8 (5.2) due a larger coupling s H _^j of the SM-like Higgs to the two pseudoscalars. 
For large m^ Q the cross section goes to zero which is consistent with the fact that this 
process is not present in the SM. We further note the thresholds in 2mt and mt + mi in 
the pseudoscalar mass, parametrized in terms of the CP-even Higgs masses according to 
eq. (5.1), become visible. These effects are not present in gg — > h^ho, since, the mass of 
the final state particles was assumed to be below these thresholds. 

5.4 Results for gg — > it (the M^-spectrum) 

In this section we present the M-mass spectrum. In the case where the h$ or po are above 
the ti-threshold {mh ,rrip Q > 2m t ) a dip-peak structure is to be expected, originating from 
the interference of the QCD-background with LW Higgs states, as described in section 4. 
This phenomenon is observed in the actual simulation as can be inferred from figure 7 
for 771/^,771^ = (125, 450) GeV but is hard to see for higher values of LW Higgs mass 
e.g. 777^,777^ = (125,800) GeV. This is because the width of the intermediate ho and po 
becomes large and the two terms in eq.(4.2) tend to cancel each other. In the latter case 
the signal to background ratio can be improved significantly in the case where a pT-cut of 
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Figure 6. Contour plot of the total cross section (in fb) for gg — > hopo (yS =7/14 TeV respec- 
tively) versus the light Higgs boson mass, rrih and heavy Higgs boson mass, for Mq = 0.5 TeV 
and Mi = M 2 = 1 TeV. 

250 GeV is applied to each top c.f. figure 7. This study could be explored further using 
the top tagger of ref. [38], since, the transverse momentum of the top quarks peak around 
300 GeV, i.e., the tops are boosted 13 . For Hq masses in the multi-TeV range one could 
employ the top tagging methods of ref. [42]. 1 Note that the two LW-states ho and po are 
necessarily close to each other in case of a low SM-like Higgs mass by virtue of the tree-level 
relation (5.1). The effect of which can be seen in figure 7 where the individual parts are 
given. We have chosen Mj r bins of 5, 15 and 30 GeV respectively for y/S = 14 TeV. A 
bin-size of 5 GeV seems unrealistic (in view of detector resolutions), whereas 15 GeV can 
be achieved and 30 GeV might very well be the reference value for early publications. A 
fundamental particle is described by its mass, spin and to some extent its interactions. 
So far we have not addressed the spin. The latter can be determined, as usual, through 
angular distributions. In [44], c.f. figure 15, the so-called Collins-Soper angle is advocated 
as promising observable. 

We would like to add that the simulations were performed with LO order QCD back- 
grounds. For an assessment of NLO corrections we refer the reader to figure 2 in [44]. 
Besides the fact that they are not too large in the low mass region the important thing is 
that the shape is very similar to LO and thus very different to a resonant structure. In 
regard to the values of the da(gg — > it)/dMi t differential cross section it should be kept in 
mind that it is not the top-pair that is observed in the detector. The efficiency of the top- 
reconstruction is estimated to be about 5% [55, 56]. The effects of the Higgs resonances for 
y/S = 7 TeV seem to small to be observed and we have relegated the corresponding plot to 
appendix A. 5 figure 14(right). In that case the gluon density is too small and qq — > g it 

13 The search strategies outlined in refs. [39-41] can, also, be applied here as well. 
14 For a review of top tagging we refer the reader to ref. [43] 
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Figure 7. Histograms of the top pair invariant mass, Mj t , for — ► for \/S = 14 TeV. (top 
left) (top right) and (bottom left) for 5/15/30 GeV-bins, respectively. A dip-peak structure is to 
observed by the interference of the QCD-background with LW Higgs states. In these figures we 
have chosen the following mass values m^ o = 450 GeV and rrip = 467 GeV which implies with 
eq. (5.1) m ho = 125 GeV. (bottom right) We plot M tt for = 800 GeV in 15 GeV bins where we 
assume Mq = M u = = 500 GeV where the signal to background ratio is significantly improved 
by f>r-cut of 250 GeV to each top. 



becomes more important. The latter being in a color octet representation, does not inter- 
fere with the LW contributions which is in a color singlet representation which leads to a 
reduction of the effect. 



6 The gg — > h h — >■ 6677 channel at the LHC 

In this section we will access the observability for double Higgs boson production in the 
LWSM 15 being the more promising than the poho-channel, for light Higgs boson masses 
in the range of ~ 120 — 130 GeV in the gg —> HqUq — > j'jbb channel. This channel is of 
particular relevance, since, searches for the SM Higgs boson at ATLAS exclude SM Higgs 
boson masses at 95% C.L. in the range 155 - 190 GeV and 295 - 450 GeV [62] and at CMS 

lj The hoho — > 7766 channel has been studied in the past in the context of the Randal-Sundrum model 
by both ATLAS [57] and CMS [58, 59], in the SM and MSSM [32, 60], and, most recently, in the context 
of a hidden sector Higgs boson [61]. 
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Figure 8. Contour plot of the total cross section (in fb) for gg — > h^ho — > blry^f (y/S = 7/14 TeV 
respectively) versus the light Higgs boson mass, mi and heavy Higgs boson mass, m~ h for Mq = 0.5 
TeV and Mi = M 2 = 1 TeV. 



exclude SM Higgs boson masses at 90% C.L. in the range 145 - 480 GeV [63] 16 . This 
suggests the SM-like Higgs boson should reside in the low mass region, i.e, mh < 145 GeV. 

Shown in figure 8 are scans at both \/~S = 7/14 TeV respectively of the cross section 
a(hoho — > 7766) over the plane of (mh ,mj io ). At = 14 TeV we choose three benchmark 
points listed in Table 1. At 7 TeV, a(hoho — > 7766) is less than or close to 1 fb throughout 
the plane of (mh , ) . This is before any sort of event selection which would reduce this 
by a factor of 10. Bearing in mind that the 7 TeV LHC is expected to accumulate about 
10 fb -1 of integrated luminosity before its upgrade to 14 TeV we do not follow 7 TeV any 
further. 

At the LHC the signal process pp — > ho — >■ /io^o — > 77&& will give rise to photons 
and jets of relatively high transverse momentum px ~ 90 GeV. In figure 9 we show 
the transverse momentum of the hardest photon and hardest jet to illustrate our point. 
Backgrounds consist of (i) di-photon plus multi-jets, (ii) single photon plus multi-jets, and 
(iii) multi-jet production. Our choice of photon isolation completely eliminates (iii) multi- 
jet production and (ii) single photon production from contention. Out of di-photon plus 
multi-jets, the dominant contributions are from the associated production of two photons 
and two heavy flavours, i.e., bottom and charm quarks. These are denoted as 77QQ where 
Q = c,b, b, c. In addition, there are backgrounds from jjQj and yyjj where j = u, d, s, g. 
Photons and jets from these backgrounds tend to be softer than those from the our signal 
process (see figure 9). 

In our simulations we model 6-tagging utilizing information in the event history of the 
Monte Carlo we are using. We label a jet a b-t&g if a partonic 6-quark of at least 5 GeV 
of transverse momentum is found in a cone of R = 0.3 around the axis of the jet. If no 

16 These bounds apply to the SM. For the LWSM we would expect, from the viewpoint of the HD- 
formalism, very similar or slightly stronger bounds. 
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6-quark is found, then we check in this order for a c-quark and r-lepton. If no heavy quark 
or lepton is found, we label the jet a light jet. Depending on which label the jet receives we 
apply the following weights: e b (E T ,rj), e mis ta g , c = 10%, e mis ta g ,r = 5%, and e mistagJ - = 0.5% 
[55, 64] 17 , which is reflected in the results in table 2. 

For the computation of the backgrounds we have applied several parton-level cuts to 
regulate any soft or collinear divergences. We require two fc^-jets with D = 0.7 and 

p\ > 20 GeV , pjj, > 20 GeV , (6.1) 
\rp\ < 2.5 , \n j \ < 2.5 , R-yj > 0.3 , # 77 > 0.3 . 

For the signal process, pp — > ho — > hoho — > 7766, we have not applied any parton-level cuts 
as there are no soft or collinear divergences. 

We simulate events at the LHC using the Monte Carlo program Sherpa 1.3.0 [65- 
68] . We have implemented the LWSM into Sherpa and have subsequently generated matrix 
elements for pp — > ho — > hoho — > 7766 using Amegic++ [69]. The matrix elements for the 
background processes have been generated using Comix [70]. All events generated include 
hadronization and shower effects. The parton shower is a Catani-Seymour subtraction 
based shower which is performed by module CSSH0WER++. Hadronization is performed by 
the module AHADIC++. Additionally, the effects of soft QED radiation off hadron and tau 
decays has been simulated using the module PH0T0NS++. 

In order to analyze events we have written an analysis plugin for Rivet 1.3.0 [71]. 
Fast j et 2.4.2 has been used to perform the clustering of final state particles into jets [72] . 
We have implemented the following selection criteria in our analysis: 

Cut 1: — Photon isolation: i) px > 20 GeV ii) pseudo-rapidity range of —2.5 < ?? 7 < 2.5 
are isolated photons if iii) J2r>r k ^T(k) < O.lpj, is satisfied where Rjk = 
\J (0 7 — <pk) 2 + (v-y ~ %) 2 an d R = 0.3. Here k can be at the particle-level 
either hadrons or photons with |?7 7 | > 2.5 or pZ < 20 GeV. 

— Exactly two isolated photons are required. 

— The hardest isolated photon is required to have a minimal transverse momentum 
of 40 GeV and i? 77 > 0.3. 

Cut 2: Exactly two fey-jets with D = 0.7 in the pseudo-rapidity range of —2.5 < rjj < 2.5 
with minimal transverse momentum, 30 GeV, are required. 

Cut 3: At least one 6-tagged jet. 

Cut 4: The di-photon invariant mass M 77 is required to be in the mass window, |M 77 — 
m ho \<2 GeV. 

Cut 5: The dijet invariant mass M&j is required to be in the mass window, |M&j — m^J < 
20 GeV. 

17 The expression for tb(ET,ri) is equal to the product of functions &e t and b v . These functions are 
explicitly shown in ref. [64]. 
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Figure 9. Shown (in arbitrary units) are the distributions for the signal process h h — »• 7766 (red) 
and one of the backgrounds, 7766 (blue), in transverse momentum of the hardest jet p^ (left) and 
hardest photon p J 1 (right). 
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Figure 10. Shown is the distribution in the invariant mass of two jets and two photons, M(,j 77 , in 
8 GeV bins for 30 fb _1 of integrated luminosity at the </S = U TeV LHC. 
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445 
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(c) 


130 
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50 



Table 1. Shown in this table are the light Higgs boson mass parameters rrih, the LW Higgs boson 
mass parameters, mi, and the mass window parameters Sm^ for benchmark points (a),(b), and 
(c). 



Cut 6: The invariant mass M^-p, is required to be in the mass window, |M^ 77 — \ < Sm^. 
Values of our choice of Sm^ o for each benchmark point are shown in Table 1. 

Table 2 displays the efficiencies and cross sections for the backgrounds before and after 
selection cuts have been applied. Efficiencies and cross sections for the signal process are 
shown in Table 3. In figure 10 we show for 30 fb _1 of integrated luminosity in 8 GeV 
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Table 2. Table of cross sections (in pb) for benchmarks (a),(b), and (c) before selection cuts (o" gen ) 
and with selection cuts (a e g) for the backgrounds QQ77, Q.7'77, and jj'77 where Q — c,b,c,b and 
j = u, u, d, d, s, s, g for V S — 14 TeV. Efficiencies (cuts 1-6) are relative where e tot is the cumulative 
efficiency. Cuts 1-3 are reproduced only once as they are the same for all three benchmarks. 
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Table 3. Cross sections (in fb) before selection and after selection for benchmarks (a) m^ = 
120 GcV, m ho = 300 GcV, (b) m ho = 130 GcV, m~ ha = 445 GcV, and (c) m ha = 130 GcV, 
m~ h = 550 GeV. Efficiencies (cuts 1-6) are relative where e to t is the cumulative efficiency. 
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Table 4. Cross sections (in fb) for h$Z — > 7766 before selection and after selection for benchmark 
(a) m,h = 120 GeV, mj io = 300 GeV. Efficiencies (cuts 1-6) are relative where etot is the cumulative 
efficiency. 
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Figure 11. Shown is the significance S/y/B + S plotted against luminosity for benchmarks (a) 
(left) and benchmarks (b) in blue and (c) in red (right). The upper and lower horizontal lines mark 
observation significances of 3ct and 5c The vertical lines represent 10 events. 



bins the invariant mass of the bjjj system for the signal scenario (a) and the sum of all 
backgrounds before cut 6 has been applied. For benchmark (a) we can expect to establish 
a 5<r-discovery with as little as 20 fb -1 . For benchmarks (b) and (c) outlook is not so 
optimistic. For scenario (b) we expect to reach 5a at 700 fb -1 and for scenario (c) we 
would need 3000 fb -1 of integrated luminosity. The primary reason for the reduced cross 
sections for scenarios (b) and (c) is that the dominant decay mode for the heavy LW Higgs 
ho is ho — > tt with Br^ Q ~ 95%. 

To this end we would like to mention that for benchmark (a) there is a background 
from ZHq production 18 . Efficiencies and cross sections are shown in table 4. It is worth 
mentioning that our analysis can be adapted for this case be changing our mass recon- 
struction hypothesis slightly. Instead of requiring the invariant mass My to be in mass 
window around the ho, we would instead, stipulate that in be in a mass window around 

18 Note that benchmarks (b) and (c) this channel is dominated by top pairs. 
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the Z boson. Additionally, the invariant mass M 77 bj should reconstruct the pq. 
7 Conclusions 

In this paper we have investigated the possibility of a light LW Higgs sector. As mentioned 
in the introduction SM-like Higgs sectors, such as the one of the LWSM, are not yet very 
well constrained as the the Higgs enters one-loop correction only logarithmically for larger 
masses and couples only very weakly to leptons obscuring the clean di-lepton detection 
channel. In practice this means that although the LW gauge bosons and the LW fermions 
are constrained to lay in the few-TeV range the Higgs sector could be very low. In view 
of indirect (EWPO) and direct (collider) constraints we have assumed the SM-like Higgs 
boson to be below then 150 GeV-value. 

We have investigated such a possibility by looking at the cross sections gg — > hoho and 
99 ~~ ^ Po^o c.f. figures 4,6 and the spectrum of gg — > it figure 7. Whereas the gg — > h$h$ 
channel is outside reach at the LHC in the SM, it is enhanced in the LWSM in the case 
where the LW-like Higgs is twice as heavy as the SM-like Higgs {m} lQ > 2ra/, ) and can 
decay at resonance through gg — > ^o^o shown in figure 1(a). The pseudoscalar 

gg — )■ pQ — > poho subprocess is close but not at resonance and turns out to be large as 
compared to SM Higgs channel but much smaller than the case discussed above as can 
be inferred from figure 6 vs 4. In our signal analysis we have therefore focused on the 
latter through gg — > hoho — > 6677 and from table 3 we see that the benchmark points (a) 
to (c) (m ho ,m~ ho ) = {(120, 300), (130, 445), (130, 550)} GeV reach 10 events for integrated 
luminosities of {8.5, 107, 595} fb -1 and the 5<r-discovery for {20, 700, 3000} fb _1 as can be 
seen from figure 11. In regard to these numbers we would like to add that the LHC is 
expected to collect 335fb _1 at 14TeV from 2012 to 2020 before the upgrade to the Super 
LHC where 1500 fb _1 is the reference number for 2025. 

The Higgs pair production cross section decreases rapidly for a ho with a mass above 
the top pair production threshold of 2mt- In this region the intermediate states ho and 
Po decay mostly into top pairs as this is the dominant decay mode, c.f. figure 15(right). 
In light of this it seems natural to investigate top pair production within the LWSM. It 
is found though that the dip-peak or in general the visibility of the resonance is diluted 
when the width is large which happens when the intermediate states can decay into top 
pairs c.f. figure 7. In the latter case the signal to background ratio can be significantly 
improved by applying px-cut of 250 GeV is applied to each top quark. An example is given 
in figure 7(bottom-right) for m^m^ = (125, 800) GeV. Further suggestions on how to 
improve the signal are given in section 6. 

Moreover, in this work we have also clarified a few things in the LWSM itself such as the 
tree-level sum rules in appendix C.l, how to reduce hyperbolic diagonalizations to standard 
methods in appendix C and the issue of spurious versus CP-violating phases in the LW 
generation Yukawa matrix in appendix C.2. Moreover we have computed box diagrams 
with two vector (gluon) and pseudo/scalar (Higgs) flavour-changing vertex analytically, 



- 21 - 



extending the results from the SM [37] and MSSM [33]. 19 The results are presented in 
appendix A. 2. 
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A Results and definitions for gg — > h h /h pQ process 

In this appendix all masses correspond to the physical masses and for the sake of notational 
simplicity we shall use the notation: 

m x>p h ya -^m x (A.l) 
for all the masses. We shall retain the subscript phys for the Yukawa matrices. 

A.l Triangle graph 

The triangle graph in the SM is given by 20 : 

^^^^ Q 

A \sM(gg ->• h h ) = 2~r^ — ^- F i/z(Pq) > Px = 4m^ phys /s (A. 2) 

S 111, j^- ~\~ ITTifJL 

where 

F 1/2 (x) = -2x(l + (1 - x)f(x)) (A.3) 

and 

{Arcsin 2 (l/-y/x) x>l 
1 A— " ( A - 4 ) 

4(MS$=l)-*0 2 x < 1 

c.f. [33] for example 21 

A. 1.1 gg —7- ho/ho — > HqHq triangles 

Since the Higgs sector (2.14) does not contribute to the loop, the LW-contribution can be 
obtained from the SM with modification of the vertices and taking into account mixing 
factors. The coupling of the Higgs to the triangle itself is modified by mixing factors in 
eq. (2.8) s H _^ and s H _^ = —s H _fj for the standard and the LW Higgs boson respectively. 
The triple Higgs boson vertices /ijj and h^h^ are modified in the same way multiplying in 



19 Flavour-changing vertices were computed in the MSSM in the squark sector [73] whereas here the top 
fermions are considered. 

20 This notation agrees with [37] as follows: a^ 2 = gaugel (2) (triangle). 

21 The function f(x) relates to the Passarino-Veltman function as follows: 2m 2 . /s (2 + (4m 2 — 
S )C (0,s,0,m 2 ,m 2 ,m 2 )) = p x {l + (1 - &,)/(#,))• 
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(a) 



(b) 




ho 



Pa 



(c) 



(d) 



Figure 12. (a-c) Triangle graphs for q = (t,i,T,b,b, B) and (d) one out of six box graphs for 
QhQj = (t,t,f,b,b,B). 

addition a factor of s 2 H ^. Furthermore Xv 2 = 2m| Q phys /(l + r h )> exceptionally insisting 
on the subscript phys, according to eq. (2.7) which leads to: 

1 1 



A (99 -> fio/ioj = 



1 + r 



/in 



s — m 



i +im ho T ho s-m 2 r-im~ h T 



I'd 



s F 1/2 (A.5) 



with 



,phys)33 j-i /o \ 

*l/2 = *l/7.\Pt) — H/2\Pi) — *l/2\Pt>) ■ 



m t 



(A.6) 



The process gg — > hopo consist of triangles and boxes shown figure 12. The triangle 
contributions can be broken down into contributions originating from: 

1. s-channel po exchange shown in figure 12(a), 

2. s-channel Zq exchange shown in figure 12(b), and 

3. s-channel Zq exchange shown in figure 12(c). 

We denote the contribution of all triangle diagrams by 

A$(gg -> h po) = A^° + A^' zt + A^' Zb , (A.7) 
where the amplitudes are further defined in the next subsection. 
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A.1.2 gg ->■ p ->• /i Po triangles 



S H-H ( 2m /i 



i > \ 2 " S ^1/2 

1 + rr \s — m 



(Ai 



where for P 1/2 = F 1/2 (F 1/2 (p x ) P 1/2 (^)) with P 1/2 (^) = AJ(/3 j;) in accordance with 
[33]. 

A. 1.3 g<? — > Zq/Zq — )• hopo triangles 



4 



A,z<? _ .eu 2 s^(cosh6'^ + sinh^) 2 



(mp Q -m fto ) x 



cos 9w sin 

3=1 



s — m| + imz^z 



s — m 2 ^ — im^Y z 




where the function / is defined in (A. 4) and = sinh(0/j) in accord with our notation 
in eq. (2.8). Note this is due to the fact that prior to diagonalization only the hopoZ- 
coupling but not the hofio Z-coupling is present. The couplings of quarks to gauges bosons 
are parametrized as follows: 



C = E (H9 Z L fI {2+^¥ L + ¥ R g Z Jf{Z+^¥ R 



f=t,b 



phys 



(A.10) 



The superscript "phys" indicates that all fields and couplings are understood to the physical 
ones. The physical couplings g^/J^y S are obtained from the expressions in Eqs. (A. 11) to 
(A.14) as 

9 L, phys = S l 9lSl , fffi.phys = Sr9rSr > j 

where X stands for Zti or Zbb respectively. 



zt l e (cosh 9z + sinh 6z ) 

UCqitSqn 



/-4(1 - cl) 



\ 



4(1 - 4) 
A ~ 2 



44 + 1/ 



(All) 



9 Z L tl 



e(cosh 6z + svcAxOz) 



Ac?,, - 1 





4(1 - 4 





(A.12) 



Zbb 

9r 



e(cosh#£ + sinhflz) 

UCqijSqn 



/2(1 -a 





\ 

-2(1-4) 

o 24 + 1/ 



(A.13) 



9 Z L bl 



e(cosh0z + smhdz) 

6C7/; Sin 



(~2cl - 1 










-2(1 - cl 







24 + y 



(A.14) 
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A. 2 Boxes for gg — > h^ho and gg — > HqPq 

For definiteness we shall give one graph, the one indicated in figure 1 (right): 

(a )? 5 )(m i ,m j )(Po)^ + ({a 2 )? 5 )(m i ,m j )(P 2 )^ |fi gure l( r ight) + X v = 
f dH 

(^ 2 i) J 7^4 tr [^ ,S 'wi( Z +Pl)lv S mi(l +Pl + P2)lS mj (l + Pi +^2+^3)75^(0] i 

for vertices 1 and 75. The term stands for are structures vanishing when contracted 
with the according polarization vectors. As stated in the main text in this notation only 
(^0,2)11 do contribute in the SM, since there are no fundamental pseudoscalars, and are 
related to the results in [37] as: (00,2)11 = gaugel(2)(6ox). 

In the following we shall present our results for the box graphs. The analytic computa- 
tions have been performed with the aid of FeynCalc [74] . We are not aware of them being 
published elsewhere for the case where the flavour can change between the Higgs vertices. 
The gluon momenta are p± and p 2 whereas the Higgs pair momenta are p 3 and are p 4 . We 
use the convention where all momenta are incoming, i.e. p\ + p 2 = ~P3 — p 4 ■ The result 
is given in terms of the Mandelstam variables 

s = {pi+ p 2 f , t = { Pl + p 3 ) 2 , u = Ox + p 4 ) 2 (A.15) 

and further shorthands 

Ti = t-m 2 , Ui = u-m} (A.16) 

for i = 3,4. 
(ao)?i(m,M) 

= -\4s + 8M 2 sC l2 + 2s((m + M){2M 2 (m + M) - Ms) - M 2 {t + u))(D 123 + D l32 + D 213 ) 
s I 

+ (m 2 3 + m\- 2(m + M) 2 ) [T 3 C l3 + T 4 C 24 + U 3 C 23 + U 4 C l4 - (tu - m 2 3 m 2 4 + s(m 2 - M 2 ))D 132 ] 
+ {m o M}| 

(a )g(m,M) 

= ^y-{- 2s(mMs + M 2 {m 2 3 - ml))(D 123 + D 132 + D 2l3 ) 

+ {ml - m\) [T 3 C l3 + T 4 C 24 + U 3 C 23 + U 4 C 14 - (tu - m 2 3 m\ + s(m 2 - M 2 ))D 132 ] 
+ {m <H> M}| 
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(00)55(771, M) = -(aoM(m,-M) = -(a )^(-m, M) 
(a 2 )°(m,M) 

= 1 9 9 ( (t 2 + u 2 - (4m 2 + 4mM) (t + u) + 4(m - M) (m + M) 3 + 2m2m 2 4 )sC 12 

tu — mim\ y 

+ (m\m\ + t 2 - 2t(m + M) 2 )(T 3 C 13 + T 4 C 24 - st^is) 
+ (m\m\ + u 2 - 2u{m + M) 2 ){U 3 C 23 + U 4 C 14 - suD 123 ) 

- (t 2 + u 2 - 2mlml)(t + u- 2(m + M) 2 )C 34 

- (i + u - 2(m + M) 2 )((tu - mlm\){m 2 + M 2 ) + s(m 2 - M 2 ) 2 )(D 123 + D 132 + D 213 )} 

+ (M 2 - m 2 )(2(m + M) 2 (u(2s + 1) - m|mf) - m 2 m 2 (s - i - u) - tu(m| + ml) - 2sn 2 )sZ?i 2 3 
+ (M 2 - m 2 )(2(m + M) 2 (t(2s + u) - m|m|) - m 2 m 2 (s - t - u) - tu(mj + m|) - 2st 2 )sD 213 
+ {m o M} 

(o 2 )°(m,M) 

= ~\ J (2(M 2 -m 2 ){u-t) -t 2 + u 2 )sC 12 

tu — mirn\ V 

+ (m 2 m 2 4 -t 2 )(T 3 C 13 + T 4 C 24 -stD 213 ) 

+ (m 2 3 m 2 4 -u 2 )(U 3 C 23 + U 4 C 14 - suD 12 s) 
+ ((t + u) 2 - 4m|mi)(t - u)C 34 

+ (t - u)((tu - m 2 m 2 )(m 2 + M 2 ) + s(m 2 - M 2 ) 2 )(D 123 + D132 + #213)} 
+i (M 2 - m 2 )((s -t + u)(tu- mlm\) + 2su(u - t)))sDi 23 
-i (M 2 - m 2 )((s -u + t)(tu- mlmf) + 2st(t - u)))sD 213 
+ {m o M} 

(a 2 )° (m,M) = -(a 2 )? 1 (m,-M) = -(a 2 )°(-m, M) (A.17) 

We would like to add three comment concerning symmetries in the amplitudes. First the 
relation, 

(ao )2 )55K M ) = "(00,2)11 K -M) = -(a , 2 )? 1 (-m, M) (A.18) 
follows from commuting the 75 from one pseudoscalar vertex to the other one. It is easy 
to see that doing this is equivalent to an overall factor of —1 and changing all the masses 
in the nominators where the 75 passed from say M — > — M. This in turn is equivalent to 
eq. (A.18). Second, the amplitudes (ao, 2 ) I [ ] 5 (m, M) can be obtained from (ao, 2 )5 1 (m, M) by 
interchanging p 3 and p 4 which results in: 

P3^P4 m 3 ^ m 4 ,u ^t ,C 13 ^ C 14 ,C 23 ^ C 24 ,D 123 ^ D 213 (A.19) 

Thirdly the aP are manifestly symmetric under interchange of t and u. We note that the 
matrix element without polzarization vectors contracted is symmetric under interchange 
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(pi,/i) <-> (p2,v) which results in t -<->• u. Thus {a) u P^ v is symmetric and since Po, P2, 
Pq (P2) are even (odd) respectively the same property holds for («o)n5/5i) > ( a o>2)(n/55) 
((°2)(i5/5i)) as can be seen from the formulae above. 

A. 3 Tensor structures 

The tensor structure for the parity-even case Po, P2 are given in [37]: 

S z = : = <r - 



O1P2) 

S = 2 ■ P» v = g^ u + P ^lP2 _ 2 (P2P3)PlP3 _ 2 (PlP3)P2P3 + V3V3 

2 Pt(^iP2) Pt(piP2) Pt{p\P2) Pt 

whereas the one for the parity-odd case [33] are: 

1 



{P1P2)P 2 



where p 2 ^ = 2(pi£>3)(p2£>3)/(piP2) — p| and the projectors {Po, Pq, P2, P2] are normalized 
as follows: 

Pi e {P , Po,P2,P2} s.t. P i P j = 25 ij . (A.20) 

Note that there are two more structures with the properties of Po an d on more with the 
property of P2. This is of no relevance as we have performed the computation by contracting 
with helicity vectors. The basis that we have chosen is p\ = {p, 0,0, p), P2 = (p, 0, 0, — p), 
e(pi,±) = e(p2,=F) = l/\/2(0, -l,=Fi,0) p 3 = (^mj + q 2 , 0, g sin(0), q cos(0)) and p 4 = 
(y / rn|~+~g^, 0, — gsin(0), — qcos(#)) where g is determined through energy conservation 2p = 
\/m\ + q 2 + -y/ m| + <? 2 . 

A. 4 Passarino-Veltman functions 

To present our results we use the standard Passarino-Veltman functions [75]: 
Cy(mi,m 2 ,m 3 ) = 

d 4 A; 1 /A . 

(A.21) 



itt 2 (k 2 — m\){{k + pi) 2 — mfydk + pi + Pj) 2 — m|) 
Ajfc(mi,m 2 ,m 3 ,m4) = 

d 4 fc 1 



i7r 2 (A; 2 - m\){{k + pi) 2 - ml)((k + pi +Pj) 2 - + +pj + j>fc) 2 - m 2 ) 

and introduce the following abbreviations 

C12 = C 12 (M, M, M) C 13 = C 13 (M, M, m) 

C u = C lA {M, M, m) C 23 = C 23 {M, M, m) 

C 2 4 = C 2i {M, M, m) C7 34 = C U (M, M, m) 

D 123 = D 123 (M, M, M, m) D 132 = D 132 (M, M, m, m) 

£213 = D 213 (M, M, M, m) . (A.22) 
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The loss of information in the exact mass dependence of the C and D functions has to be 
taken into account when symmetrizing in m and M in formulae Eqs (A. 17). 

A. 5 Additional plots 





Figure 13. Contours plots of the ratio Xh = log ( ^r 2- ) (right) of the ho and the ratio x% 
log ( of the ho. 




Figure 14. (left) Width as a function of mass, m ho , for m, iQ = 120 GeV, M 2 = Mi = 1 TeV 
for different values of the fermion mass scale, (right) Histogram for gg -> for = 7 TeV with 
5 GeV-bins. 
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Figure 15. (left, right) Branching ratios Br/^ and Br^ Q as a function of the masses m^ and m^ 
and fixed m~ ho =120 GeV and m ho = 450 GeV respectively for M 2 = M t = 1 TeV and Mq = 500 
GeV. 




i i 1 ' 1 1 1 ' 1 1 i i 1 1 1 ' 1 1 1 1 

200 300 400 500 600 700 800 900 1000 200 300 400 500 600 700 800 900 1000 
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Figure 16. The cross section of pp — > h n h via gluon fusion at the LHC for yfs — 7/14 TeV 
respectively versus the mass of the ho, m^ o , for m/i = 120 GeV. We note that the fermion mass 
scale Mq = M u = Md has very little influence on the results as emphasized in section 5. Note that 
for large the SM model value is approached by virtue of decoupling of the LW Higgs. 



B Results for gg — > h /h /p — >■ tt 

The amplitudes for the processes can directly be obtained from the ones from the double 
Higgs pair production in the previous section by suitable replacements. Prom the amplitude 
gg — > ho — > hoho in eq. (A. 5), using eq.(2.14) and the definition of A in the Higgs potential 
chosen in section 2.1 one obtains: 

Ao(gg -> h (h ) -> tit) = s 2 & {gl h ) u ( 3— —)s F 1/2 [it] 

Furthermore, from the gg — > po — > poho amplitude in eq. (A. 8) one obtains: 

A^°(gg ^po^ it) = ( ~ 2z(g P^ )n s p ) M (B.l) 
Note, in both cases, we have not evaluated the spinors t, i. 
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C Diagonalization of Mass Matrices 



Here we shall describe a method for performing the hyperbolic diagonalization 

M t , physV3 = S ] R M t mS L (C.l) 

using similarity transformations for which standard tools, e.g. Diag 1.3 [76], can be used, 
based on the observation that: 

{s R/L mr l = s RJL m (c.2) 

The latter relation is easily verified from eq. (2.13) 

Here we will describe a procedure of obtaining Si and S R numerically using routines 
provided. From there it is straightforward to verify that: First, we recognize that 

diag(m t 2 phys ,m? phys ,m| phys ) = Mt, ph mMl pi ji3 

= A R { m M tm M\)A- R x = Al^mIvsMJAI 1 . (C.3) 

with A R = S R r] 3 and A L = 
C.l Mass sum rules 

In this section we would like to point out some tree- level sum rules for matrices. When 
the matrices are diagonalized by hyperbolic rotations the trace remains an invariant. To 
be more precise suppose we had a matrix that is diagonalized as follows 

M physV = S^Mr ] S , (C.4) 

with 

S^S = r) , M phys = diag(m* phys , mf phys , ...) , (C.5) 

then 

trlMphys] = tr [-M] • (C6) 

The correctness of (C.6) can be immediately verified using the properties above. The 
diagonalization can be interpreted as a symmetry transformation where r/ plays the role 
of the metric. Thus the statement eq. (C.6) is nothing but the fact that the trace of the 
(2,0)-tensor (Mrj) a p is an invariant; M Q a = tr[M]. Thus one can deduce sum rules for the 
masses. Applied to the CP-even Higgs sector the RHS follows from writing (2.4) in matrix 
form, c.f. [26] and the LHS is given by definition 

m Lphys + m i = M 2 H = (m| 0jphys ) . (C.7) 

The correctness is readily verified from eq. (2.4). Note, with the peculiar fact that at tree- 
level m~ o phys = Mft equation (5.1) follows. This technique applies to the entire bosonic 
sector. For the neutral gauge bosons one gets 

< phys + ™ W + ™|, phys = Ml + Ml (C.8) 



- 30 - 



with M\p the mass scale of the U(l)y and SU(2)i HD gauge terms respectively. The field 
A is the LW-partner of the photon. Note the photon is not explicitly written down since 
it remains massless. eq. (C.8) is consistent with the result for M\ = M2 in appendix B of 
reference [24]. 

The fermions are slightly more complicated as they proceed via a bi-unitary hyperbolic 
diagonalization. The statement is that: 

diag(m 4 2 phys , m| phys , m\ phyg ) = tr[A^ ph% Mj )ph 77 3 ] = tx[M t mM\m] , (C.9) 

which follows immediately from the eq. (C.3). Applied to the fermions we get: 

m ?, P hys + m| phys + m 2 f phys = Ml + Ml , (CIO) 

where eq. (2.11) was invoked for A4f The correctness of this equation can be verified for 
the explicit result given in chapter 2.3.2. of reference [26] to each order in the expansion. 
In chapter 3 of reference [26] similar consideration were taken into account to show the 
absence of quadratic divergences in the top-loop in the AF formalism. 

We would like to emphasize that the trace formula (C.6) and (C.9) are general and in 
particular apply in each order of perturbation theory but the specific evaluation we have 
given in Eqs (C.7),(C.8) and (C.10) have made use of the trace at tree-level and are thus 
subject to corrections. 

C.2 Spurious phases 

Furthermore we consider it worthwhile to discuss the freedom of reparametrizing phases 
in the mass and Yukawa matrix of the LWSM. Note that the Yukawa matrix presented 
in ref. [26] contains imaginary entries and one might therefore wonder whether they are 
associated with CP-violation or whether they are unphysical/spurious phases. For fixed 
flavour there are six fermion in each LW-generation counting left and right handed field 
separately. The freedom of choosing their spurious phases is reflected in the fact that the 
matrices Al and Ar are determined by eq.(C.3) up to 

A R -> &a,g{e iR \e lR \e iR3 )A R , A L -> diag(e iL \ e iL \ e iL[i )A L (C.ll) 

a multiplicative diagonal unitary matrix. Rewriting eq.(C.l) as 

M ttPhys m = AnirjsMmUi 1 (C.12) 

we see that choosing the fermion masses to be real and positive (or negative) fixes the 
differences R{ — Li for i = 1,2,3. Writing L\ = L\,L2 = L± + A2,L 3 = L\ + A3 it is 
noticed, as usual, that only the two parameters A2 and A3 lead to a change in the entries 
of 

9t,phys'i two arbitrary phases. This freedom can be used to reparametrize the third LW- 
generation by e %Rs = e tLi = i the Yukawa matrix gt,phys m r ef. [26] to render its entries 
completely real. 

To this end we would like to note that we find that <?t, p hys is smooth in the limit 
Mq — > M u contrary to a remark made in the appendix of ref. [77]. Note in their explicit 
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formula these authors present an expansion in 1/(M U — Mq) which cannot be compared 
with the expansion in 1/M U for M u = Mq presented in ref. [26] of as the former is singular 
in the degenerate limit. The fact that their expansion does not have imaginary parts can 
be explained by the freedom of phase reparametrization discussed above. 
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